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Abstract 

We represent the Polyakov loop correlator as a spectral sum of cor- 
relators of eigenvectors of the lattice Dirac operator. This spectral rep- 
resentation is studied numerically using quenched SU(3) configurations 
below and above the deconfinement temperature. We analyze whether 
the individual Dirac eigenvector correlators differ in the confined and de- 
confined phases. The decay properties of the normalized Dirac eigenvector 
correlators turn out to be essentially identical in the two phases, but the 
amplitudes change. This change of the amplitudes shifts the relative con- 
tributions of the individual Dirac eigenvector correlators and is the driving 
mechanism for the transition from the confining static potential into the 
deconfining one. 



Introductory remarks 



Confinement of quarks is one of the key features of QCD. Of particular interest is 
the fact that at a critical temperature QCD undergoes a transition into a decon- 
fined high temperature phase. Understanding the transition and characterizing 
QCD in the deconfined phase are widely discussed problems. 

It is obvious, that confinement and deconfinement, respectively, must be 
manifest in the Dirac operator and its inverse, the quark propagator, since the 
quarks need to know whether they are confined or not. In a series of recent 
papers [lj - |7j the connection of quantities related to confinement, e.g., the 
expectation value of the Polyakov loop, to spectral sums of the Dirac operator 
was studied. Signatures of confinement in the Dirac spectrum and their change 
at the QCD transition were identified. 

In this article we go beyond results for only the Dirac eigenvalues and study a 
spectral representation of the static quark anti-quark potential. Since this quan- 
tity depends on the relative spatial position of the static sources, its spectral 
representation contains also the eigenvectors of the Dirac operator. In particular 
the shape of the potential must reflect itself in spatial correlators of Dirac eigen- 
vectors. Very little is known about such correlators of Dirac eigenvectors and 
we study these relations numerically using quenched lattice QCD configurations. 
We explore how the spectral signatures of Dirac eigenvalues and eigenvectors 
change at the deconfinement transition. 

We remark at this point that very recently Synatschke, Wipf and Langfeld 
have studied correlators of spectral sums with a different, IR dominated spectral 
function [7]. Using the boundary condition technique of [lj they project to loops 
with the same transformation properties as the Polyakov loops which we use 
here. The continuum arguments given in [7J and the numerical study with SU(2) 
gauge configurations provide interesting complementary insight to the analysis 
presented here. 



Dirac operator and Polyakov loops 

In this section we repeat and refine the arguments [1] which lead to the con- 
nection of the Polyakov loop and its correlators to spectral sums for the Dirac 
operator. 

The Polyakov loop at spatial position x is given by 
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U4(x,t) is the temporal component (= 4-component) of the gauge variables 
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Ufj,(x,t) £ SU(n) and tr c is the color trace. The rank n of the gauge group 
in the fundamental representation appears in the overall normalization of the 
Polyakov loop. The product runs over all N temporal links at a given spatial 
point x. We work on a lattice of size L 3 x N and for the gauge field use 
boundary conditions which are periodic in all four directions. We are interested 
in correlators of the Polyakov loop 

(L(x)L*(y)) ~ Ae- NV U*-M, (2) 

where * denotes complex conjugation and V(\x— y\) is the static quark-antiquark 
potential. 

The staggered Dirac operator at vanishing quark mass is given by (we set 
the lattice spacing to 1) 

1 4 

D (x, y) = r ^( x ) u ^ x ) ~ U v( x ~ A) f *x-A,i/ , (3) 

z v=l 

where x and y are integer valued 4-vectors labeling the lattice sites and rj^x) = 
(— l) Xl+ ■■■ js the staggered sign function. 

The terms of the staggered Dirac operator Q connect nearest neighbors. 
When powers of D are considered, these terms combine to chains of hops on the 
lattice. Along these chains products of the link variables U^(x) are collected. 
Taking the m-th power will give rise to chains of length m. We now consider the 
iV-th power of D, where N is the temporal extent of our lattice. Thus we will 
encounter chains with a length of N. Furthermore we set the two space-time 
arguments of D to the same value, y = x, such that we pick up only closed 
chains, i.e., loops starting and ending at x. Among these are the loops where 
only hops in time direction occur such that they close around compact time. 
We obtain 

(_l)N(x 1 +x 2 +x 3 ) N 

tr c [D N (x, t\x, t)] = trivial loops + i '— „ tr c U 4 (x, s) 

1 8=1 

(_l)N(x 1 +x 2 +x 3 ) N-l 

+ ^ ' 9jV tr c II M^N-S)^ 

Z s=0 

= trivial loops + ^raL(x) + ^nL*(x) . (4) 

The factor (— l) 7 ^ 1 +^2+^3) comes from the products of the staggered sign 
factors 7/4(2;). Here we use even N, such that this factor is equal to +1. We 
stress that the forward and backward running Polyakov loops are the only ones 
that wind non-trivially around compact time. 
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We now explore the fact that the Polyakov loops respond differently to a 
change of the boundary conditions compared to other, trivial (= non-winding) 
loops. We can change the temporal boundary condition of the Dirac operator 
by multiplying all temporal link variables at t = N with some phase factor 
z g C, |z| = 1: U±{x,N) -> zU±(x,N) Vf. We evaluate the left-hand 
side of (TjJ using the Dirac operator in the transformed field, and denote the 
corresponding Dirac operator as D r , where r labels the boundary conditions 
we use. Actually we linearly combine three different boundary conditions z r , 
r = 1, 2, 3 to obtain: 

2 N 3 3 

— c r tr c [D^(x, t\x, t)] = L(x) x^CfZy 

lb -. 

r=l r=l 

3 3 

+ L*(x) x 2_. Cr z* + trivial loops x YJ c r . (5) 

r=l r=l 

The complex coefficients c r are determined by the requirement that the factors 
in front of the complex conjugate loop L{x)* and in front of the trivial loops 
vanish, while the factor in front of the Polyakov loop L(x) is required to be 
equal to 1. Here we use the following values for the z r and the c r , 

Zx = 1 , z 2 = e i2n/3 , z 3 = e~ i2n/3 and c r = z*/S , (6) 

which fulfill the above requirements. We conclude 

2 N N 3 

L ( f ) = ^ E E < ^ [D?(x, t\x, t)\ . (7) 

t=l r=l 

We have used the fact that the left-hand side is independent of t and averaged 
the right-hand side over t. 

Spectral sum for the Polyakov loop correlator 

We now express the N-th power of the Dirac operator as a spectral sum. Let 
v^ k \x,t) be the eigenvector of D with eigenvalue \( k \ Then 

D N (x,t\x,t) c , c , = Y,(^ k) ) N v W (x,t) c v^(x,ty c , , (8) 

k 

where we have made explicit also the color indices c, d . The index k runs over 
all eigenvalues of the Dirac operator. Inserting the spectral sum into (T7]) we find 

2^ 3 N 
= rEE< fl P^(x), (9) 

k r=l 
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-i N n 2 

p { r ] ^) = jtY,Y\$ k) &M > ( 10 ) 

t=l c=l 

where again r labels the three different boundary conditions we use and the 
index k runs over all eigenvalues. 

The spectral sum may be simplified further, using the fact that the eigenval- 
ues of the staggered Dirac operator are purely imaginary and come in complex 
conjugate pairs. The eigenvectors corresponding to these pairs are related to 
each other by multiplication with the staggered sign (— i^i+^+a^+o;^ which 
leaves the densities (llOp invariant. Thus in the spectral formula one needs 
to sum over only the eigenvalues with positive imaginary part (we indicate this 
with the superscript + attached to the summation symbol) and an extra factor 
of 2 appears. Thus we obtain for the Polyakov loop correlator 

A N + l 4- 3 N 
L(x) L*(y) = -g^-E E 4*. (AW) P?\2) P { JHy)- (11) 

k,l r,s=l 

We stress that the result (fTTj) is an exact formula, in particular it holds for 
an individual gauge field configuration and not only in the ensemble average. 
Furthermore, all terms involved in our spectral representation are manifestly 
gauge invariant quantities. 

Let us finally discuss the correlator after the vacuum expectation value (..) 
has been evaluated. In this case one can use the fact that the expectation value 
is invariant under spatial reflections, i.e., we can interchange x and y. When 
using the choice ([6]) for the phases at the boundary one obtains a particularly 
simple form of the correlator 



, aN+1 



W^A^V'W^)) (12) 



N 



The final result is a sum of 6 terms, where the 3 correlators of the densities 
with equal boundary condition come with a plus sign, while the correlators with 
mixed boundary conditions are subtracted. 



Truncated spectral sums 

The formulas (jlll) and f 1 1 2 p constitute the announced spectral representation of 
the Polyakov loop correlator and thus, via Eq. (J2j), the spectral decomposition of 
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the static quark-antiquark potential. The right-hand side is a double sum over 
all Dirac eigenvalues and we can indeed decompose the correlator into IR and 
UV parts by taking into account only parts of the spectrum. We now analyze 
numerically how the individual terms in pip and (1121) give rise to the static 
quark potential and how they see the QCD phase transition. 

The numerical analysis is done on quenched SU(3) gauge configurations 
generated with the Liischer-Weisz action [8] on various L 3 x N lattices with 
L = 10 and 12 and N = 4 and 6. All data we show here are for the largest 
lattice, 12 3 x6, and we have a statistics of 20 configurations both in the confined 
and the deconfined phase. These results were cross-checked on smaller lattices 
where larger statistics of 0(100) configurations were produced. The lattice 
spacing was determined in [9J using a Sommer parameter of vq = 0.5 fm. The 
gauge coupling was chosen such that we have ensembles in the confined phase 
at a temperature of T = 236 MeV and ensembles in the deconfined phase at 
T = 412 MeV. Above T c (~ 300 MeV) the theory has a non-vanishing Polyakov 
loop, which can have three possible phases. In our numerical analysis we restrict 
ourselves to configurations with (essentially) real Polyakov loop, since the results 
for ensembles with one of the two complex phases of the Polyakov loop may be 
related to the real case by a simple Z3 transformation. For the numerical analysis 
of the spectral sums (TTTT> and (fT2j) we computed complete sets of eigenvalues 
and eigenvectors of the staggered Dirac operator using ScaLAPACK routines. 
The results from the spectral sums were cross-checked to an evaluation of the 
Polyakov-loop correlators directly from the gauge links and agreement of the 
first 7 to 10 digits was found. 

The central quantity we consider in this section is the correlator (L(x)L*(y)), 
evaluated with the exact formula (TTTI) . In order to improve the statistics we av- 
erage over x and y, keeping the distance d fixed. We define 

CV) = TfsE E (mL^x + dej)), (13) 

x .7=1,2,3 

where = 1,2,3 are the unit vectors in the three spatial directions and 
for L(x)L*(x + dej) we insert the spectral representation (TTTTl . According to 
Eq. © C(d) behaves as 

C{d) ~ Ae- NV{ $ . (14) 

We now study the correlator C(d) but sum the spectral sums in (fTTj) only 
up to some largest eigenvalue X cut . Thus for small \ cut we only see the IR 
contributions and as we increase X cu t towards the largest eigenvalue X m ax, we 
gradually recover the full result. These truncated sums allow one to study how 
the full correlator builds up from the different contributions in the IR and UV. 
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Figure 1: Absolute values of the truncated spectral sums for the Polyakov 
loop correlator C(d) for our 12 3 x 6 lattices below (lhs. plot) and above T c 
(rhs.). The symbols with error bars represent the full correlator as obtained 
directly from the gauge fields. The other curves are obtained from the 
spectral sums at various values of X cu t, starting from a small cutoff for the 
curves at the bottom. 

In Fig.[T]we show the absolute values of the truncated sums for our ensembles 
on 12 3 x 6 lattices below (lhs. plot) and above T c (rhs.). The symbols represent 
the data from an evaluation of the correlators C(d) using directly the gauge 
links, and thus are the limiting values which the spectral sums approach as 
the cutoff \ cut is sent towards the largest eigenvalue. Only for those we show 
statistical errors, which are however representative also for the truncated spectral 
sums. In the confined phase, the correlators show a cosh-type of behavior in d 
which is a consequence of a non-vanishing string tension for the potential V(d) 
in ([2j) and the periodic boundary conditions for the gauge fieldsj. For T > T c 
(rhs. plot) the symbols quickly form a constant plateau as is expected in the 
deconfined phase where the string tension vanishes. 

x We remark, that the rather small lattices which we are restricted to by the need to 
compute the full Dirac spectrum lead to sizable statistical errors towards the minimum of 
the cosh. The error bars appear asymmetrically due to the logarithmic scale. 
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The curves without symbols and error bars are the results from the spectral 
sums at different values of the cutoff X cut , with the curves at the bottom of the 
plots corresponding to a low X cut . We remark that due to the different values of 
the gauge coupling which we use below and above T c , a different lattice cutoff 
results, such that for the ensemble below T c (Ihs. plot) the largest eigenvalues 
are near A max ~ 3500 MeV, while above T c (rhs.) they reach up to A max ~ 5800 
MeV. Both phases show a similar behavior of a slow buildup of the correlator 
as \ cut is increased. At about 30 % of the maximal eigenvalue (e.g., the \ cu t = 
1125 MeV curve for the Ihs., and the \ cut = 1575 MeV curve on the rhs.) 
the shape of the correlator, cosh-type and flat, respectively, is quite visible. As 
X cut is increased further, fluctuations appear which might be a consequence of 
the fact that contributions from large eigenvalues are weighted much stronger. 
Towards the UV end, the spectral sums even overshoot the final result and settle 
at the exact values (symbols in the plot) only when all eigenvalues are taken 
into account. This overshooting phenomenon was observed already in |2J where 
the expectation value of a single Polyakov loop was studied. 

Spatial structure of the eigenvector density correlators 

Having established that the spectral sums give rise to the correct Polyakov loop 
correlator and carry the information about confinement versus deconfinement 
also when not fully summed, we now come to analyzing individual parts in the 
spectral sums. In particular we discuss the role of the eigenvector densities 
pf\x) in the two formulas flTTJ and (fT2"|) . 

The correlator (fT2l) is related to the static quark potential as stated in 
Eq. ([2]). Below T c the theory is confining with a non-zero string tension a, such 
that for sufficiently large \x — y\ we expect (L(x)L*(y)) cx exp(— Na\x — y\). 
Above the critical temperature the confinement is gone, such that (L(x)L*(y)} 
is essentially constant for large \x — y\. The characteristic spatial behavior of the 
correlator (fl2|) below and above T c can only come from the eigenvector densities 
p\ k \x), since their correlators are the only terms that contain a dependence on 
the spatial coordinates x and y. Consequently the correlators of the eigenvector 
densities must already contain the spatial information which in the full correlator 
(fT2l) gives rise to the discussed spatial behavior governed by the static potential 
below and above T c . 

In order to study the spatial structure of the (weighted) correlators of the 
eigenvector densities, we analyze 

C ^( d ) = E ((^r^K^P^ix + de,)). (15) 

6L x j=l,2,3 
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A and \i are eigenvalues on the positive imaginary axis and p r (x),p s (x+dej) 
are the corresponding eigenvector densities. The indices r,s = 1,2,3 refer to 
our three boundary conditions z%, Zi,z% as stated in ([6]). 

The correlators (fl5j) are rather involved objects. Beside the dependence 
on the distance d, they have four different parameters: two labels A, [i for 
the eigenvalues and two labels r, s for the 3 different boundary conditions. 
Thus we have to consider correlators where A = fi, as well as correlators with 
A ^ /i. Furthermore our formulas (fTTj) , (fT2j) contain all possible combinations of 
boundary conditions. In particular in the formula for the ensemble average (fl2|) 
the correlators with equal boundary conditions, r = s, come with a plus sign, 
while the correlators with differing boundary conditions, r ^ s, are subtracted. 
All these possible combinations of the parameters A, (J,, r, s, might give rise to 
a different dependence of the corresponding correlators on the spatial distance 
d. In (fT2l) these correlators are combined with relative signs and only this 
combination then gives the correct dependence of the Polyakov loop correlators 
on the spatial separation d. 

The central question we want to answer in this section is whether also the 
individual eigenvector density correlators show a change in their spatial behavior 
as one crosses from the confined into the deconfined phase. A changing spatial 
spatial behavior of the eigenvector density correlators would be the simplest 
mechanism for changing the static potential from confinement to deconfine- 
ment. 

While for the behavior of Dirac eigenvectors with small eigenvalues quite 
some information is available, very little is known about Dirac eigenvectors with 
eigenvalues that are not in the deep IR. Thus we here systematically study the 
eigenvector density correlators (fl5|) for various combinations of the parameters 
A, //,r, s. Since we here are predominantly interested in the decay properties 
of these correlators, we normalize them such that they assume the value 1 at 
vanishing distance d = 0. 

In Fig. [2]we begin with the simplest case, where we consider the correlation 
Cr*r' X \d) of an eigenvector with itself, i.e., both densities p^\x) in the corre- 
lator correspond to the same eigenvalue A = fi and have the same boundary 
conditions which we set to r = s = 1, i.e., we consider periodic boundary con- 
ditions for the Dirac operator. We remark, that the other two cases where the 
boundary conditions are equal, r = s = 2 and r = s = 3, give similar results. 
The Ihs. plot shows the correlator below T c , while the rhs. is for the deconfined 
phase T > T c . On the axis from the front to the back the spatial separation 
d is plotted, while on the axis running from left to right, we plot the number 
of the eigenvalue A, where at the left we start with the IR end of the spectrum 
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Figure 2: Normalized correlator Cr,r yd) for equal eigenvalues as a function 
of d and A. The boundary conditions are chosen periodic, i.e., r = 1. The 
lhs. plot is for T < T c , while the rhs. is for the deconfined phase above T c . 

and the UV eigenvalues are at the right-hand side. Note that the eigenvalues 
are purely imaginary and each eigenvalue comes with both signs. As we have 
already discussed it is sufficient to sum over only one half of the eigenvalues, 
and in the corresponding axis of the plot we simply run through the eigenvalues 
on the positive imaginary axis. These eigenvalues were divided into bins, and 
the correlators averaged over the bins is what is shown in the 3-D mesh plot. 

For both cases, T < T c and T > T c the correlator decays with d, which 
due to the periodicity gives rise to a minimum in the center followed by another 
increase. The decay is most prominent for the eigenvalues at the IR and UV 
ends of the spectrum, while near the middle of the spectrum a rather flat trough 
is formed. A more detailed analysis of slices through the 3-D plot reveals that 
at the IR and UV ends the correlator indeed gives rise to a cosh, which shows 
that the correlator decays exponentially. Slices through the center part of the 
plot show that for the bulk of the spectrum the correlators of densities with 
medium sized eigenvalues become flat. 

However, the most remarkable finding of the plots is the fact, that the nor- 
malized density correlators below and above T c are almost indistinguishable. In 
Fig.[2]this is only a finding for a very particular correlator, namely the one where 
the eigenvalues of the two densities and their boundary conditions coincide, but 
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Figure 3: Same as Fig. [21 but now the correlator C\ 2 ' (d) which combines 
periodic conditions and the density with boundary phase 2ir/3 is shown. 

we will show, that for essentially all normalized correlators contributing in (fTT!) 
and (fl2l) the difference between the confined and deconfined phases is very 
small. 

The next correlator of densities we study is the combination where the two 
densities still correspond to the same eigenvalue, but have different boundary 
conditions. In particular we consider the combination of a periodic boundary 
condition and one with phase 2ir/3. The normalized correlator Ci% ((£) for 
this choice of parameters is shown in Fig. [3J As in the case of equal boundary 
conditions we see that a pronounced decay is visible only for the smallest and 
the largest eigenvalues, while for the eigenvalues in the center of the distribu- 
tion again a rather flat plateau is found. Also here the difference between the 
confined (Ihs. plot) and the deconfined phase (rhs.) is rather small. 

So far we have restricted ourselves to correlators where the two densities 
were for equal eigenvalues. In Fig. [4] we now show some cases of normalized 
correlators of densities Cr,«'^(d) corresponding to different eigenvalues, A ^ \i. 
We always keep the eigenvalue fj, fixed and plot the correlator as a function of 
the other eigenvalue A and the spatial separation d. For the Ihs. column /j, is 
chosen near the IR end of the spectrum, while on the rhs. fi is at the UV end. 
In the top row both densities are for periodic boundary conditions, while in the 
bottom plots periodic, combined with 2n/3 boundary conditions are used. All 
plots in Fig. [4] are for T < T c . It is obvious, that these correlators mainly give 
rise to a flat spatial behavior. Again the corresponding correlators for T > T c 
are almost identical. 

From our analysis of the spatial behavior of Cr,« ,A ^(d) we conclude, that the 
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T < T 



T < T c 




Figure 4: The normalized correlators Cr,s \d), where /i is held fixed at 
an IR value (lhs. column of plots) and at an UV value (rhs.). In the top 
row we show the correlator where both densities are for periodic boundary 
conditions, while in the bottom row, the combination of periodic with 27r/3 
boundary conditions is displayed. All plots are for T < T c . 

eigenvector density correlators (fl5|) show a different spatial decay depending on 
the chosen combination of the parameters \,fJ,,r and s. For some values we 
identify a clear exponential decay, which due to the periodicity is manifest as 
a cosh, while for other choices of the parameters the correlators are essentially 
flat. The final formula (fl2l) for the Polyakov loop correlator contains a sum 
over all parameter combinations. The individual correlators then combine their 
specific different spatial shapes such that either the exponential decay for the 
confined phase, or the flat behavior of the deconfined phase emerges in the sum 
for the complete Polyakov loop correlator. 

The most remarkable observation for the normalized eigenvector density 
correlators is, however, the fact that they are almost unchanged if the results for 
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T <T C and T > T c are compared. We quantified this observation by analyzing 
the difference of the correlators ci X s ' IJ '\d) for our ensembles below and above 
T c . We found that the relative difference is below one percent throughout. The 
largest discrepancies (about one percent) were found at the IR and UV edges 
of the distributions. 

We repeated the same analysis for the modified correlators 

C^\d) = E (pPWpPP+dSjj), (16) 

x j'=l,2,3 

where the weight factor {\n) N is omitted. Although these eigenvector density 
correlators are not directly related to the Polyakov loops, their behavior is of 
general interest: Hadronic correlation functions may be expressed with the spec- 
tral theorem and various correlators of Dirac eigenvectors will appear in such a 
representation. When crossing the phase transition, the changing hadron con- 
tent might leave a trace in the correlators of the Dirac eigenvectors. In our 
analysis of the correlators Cr^\d) we found that their general behavior is very 
similar to the one of the weighted correlators Cr^\d) as shown in Figs. 2-4. 
By analyzing as before the difference of the correlators above and below T c we 
find again, that crossing the phase transition leads to changes of the correlators 
which are only in the one percent range. 

Amplitudes of the density correlators 

In the last section we studied the spatial shape of our correlators (fT5l) and thus 
plotted them such that they were normalized to 1 at d = 0. We now focus 
on the actual amplitudes of these correlators defined as the correlator Cr,a (0) 
taken at distance d = 0. These amplitudes are important for understanding how 
the individual correlators combine in the formula (fT2jl to build up the Polyakov 
loop correlation function. 

In Fig. [5] we plot the amplitudes Cr,s'^(0) as a function of A and \i and 
compare different combinations of the boundary condition parameters r and 
s. The plots on the Ihs. are for boundary conditions r = s = 1, and on the 
rhs. the combination r = l,s = 2 is shown. The top row is for T < T c , while 
at the bottom T > T c is shown. The plots show that the amplitudes for the 
contributions of different eigenvalues all have a very similar size, and only at the 
IR end show a sharp decrease. Furthermore the two combinations of boundary 
conditions (r = s = 1 in the Ihs. plot and r = 1, s = 2 on the rhs.) give rise to 
almost indistinguishable amplitude distributions. The same is true for the other 
combinations of boundary conditions that enter in the spectral sum (jl2|) . 
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Figure 5: The amplitudes C# w (0) function of the eigenvalues A and 
fi. The plots on the lhs. are for the boundary conditions r = s = 1, while at 
the rhs. we show r = 1, s = 2. In the top row we show the amplitudes for 
T <T C , while at the bottom we have T > T c . 



If we now compare the amplitudes below T c (top plots) to their counter- 
parts above T c (bottom), we find a sizable discrepancy for some ranges of the 
eigenvalues A and fi (note that we use a logarithmic scale on the vertical axis). 
To illustrate this shift of the amplitudes more clearly we define the relative 
amplitudes 



cft M) (o) 



T<T C 



(0) 



T>T C 



^r,s 



(0) 



T<T C 



+ C&'ri (0) 



(17) 



T>T C 



The relative amplitudes Ar jS may assume values between —1 and 1 and are 
positive when the correlator amplitude is larger in the confined phase, while it 
is negative otherwise. For identical amplitudes one has A^'^ = 0. 
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In Fig. Owe show the relative amplitudes A^ { on the Ihs. and A^ £ on 
the rhs. The two different combinations of boundary conditions give rise to 
very similar behavior. Both relative amplitudes have in common that they are 
negative except when both eigenvalues A and /i are in the upper third towards 
the UV end of the spectrum. This implies that the amplitudes in the IR are 
smaller in the confined phase, while towards the UV part of the spectrum the 
amplitudes are larger in that phase. Thus in the confined phase more weight 
is on the UV end of the spectrum, while in the deconfined phase the IR part 
is enhanced. This shift of the amplitudes changes the weight of eigenvector 
density correlators with different decay properties, which we analyzed in the 
last section, and in this way switches from the confining to the deconfining 
static potential. Since the spatial decay properties of the eigenvector density 
correlators are essentially invariant, it is this relative shift of the amplitudes 
which we identify as the mechanism giving rise to the changing behavior at T c , 
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Conclusions 



In this article we have studied correlators of Polyakov loops and their spectral 
decomposition in terms of eigenvalues and eigenvectors of the Dirac operator. 
The spectral sums we discuss connect the static quark-antiquark potential to 
correlators of densities of the Dirac eigenvectors. From those eigenvector density 
correlators the static potential must build up its spatial dependence, character- 
istic for the confined and deconfined phases. The eigenvector density correlators 
were studied numerically for quenched SU(3) configurations below and above 
the transition temperature T c . 

We find that depending on the parameters (eigenvalue combination and 
boundary conditions), the individual eigenvector density correlators have differ- 
ent decay properties, ranging from a strong exponential decay to a rather flat 
behavior. The final Polyakov loop correlator is obtained as a combination of all 
these contributions and only for the full spectral sum the spatial decay properties 
of the Polyakov loop correlator emerge. 

An interesting outcome of the numerical part of our analysis is the finding 
that the spatial decay properties of the (normalized) correlators of the eigenvec- 
tor densities are essentially the same in the confined and the deconfined phases. 
The relative difference of the normalized correlators below and above T c was 
found to be less than one percent. For the amplitudes the situation is different. 
We show that they change when crossing T c , such that in the deconfined phase 
the IR part is enhanced in the spectral sums, while the UV part is reduced. This 
shift of the amplitudes leads to a different combination of spatial decay proper- 
ties of individual correlators and the static potential switches from confinement 
to deconfinement. 

Having established a particular behavior of Dirac eigenvector correlators in 
the two phases as mechanism for shifting from confinement to deconfinement, a 
more systematic analysis of general eigenvector correlators appears interesting. 
Since with the spectral theorem hadronic observables may be expressed in terms 
of Dirac eigenvector correlators, a better understanding of such correlators might 
lead to a clearer characterization of hadronic properties in the deconfined phase. 
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